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THE  EFFECT  OF  WAVE  BREAKING  ON  THE  DIRECTIONAL 
SPECTRUM  OF  WAVES  IN  WATER  OF  VARIABLE  DEPTH 
IN  THE  PRESENCE  OF  CURRENT 

PART  1 :  INTRODUCTION 

1 .  Waves  undergo  changes  as  they  propagate  from  deep  to  shallow  water  and/or 
as  they  meet  an  adverse  current.  The  factors  that  affect  these  changes 
include  input  from  local  wind,  non-linear  wave-wave  interactions  and  wave 
breaking.  In  this  study,  only  the  effect  of  the  last  factor  is 
considered . 

2 .  There  exists  a  number  of  methods  to  compute  changes  of  wave  character¬ 
istics  under  the  conditions  described  above.  One  analytical  method  is  to 
use  energy  balance  equation  (Phillips,  1980,  Huang  et  al,  1972)  but  the 
method  does  not  take  into  account  the  effect  of  wave  breaking. 

3.  In  a  previous  study  (Tung  and  Huang,  1987),  a  method  is  developed  to 
study  the  effect  of  wave  breaking  on  spectrum  of  random  waves.  The  waves 
are  assumed  to  be  long-crested,  travelling  normally  toward  a  straight 
shoreline  over  gently  varying  bottom  with  straight  contours  parallel  to 
the  shoreline,  where  they  encounter  a  horizontal  current,  steady  in  time 
and  non-uniform  in  the  horizontal  space,  where  flow  velocity  is  parallel 
to  the  direction  of  the  waves  and  is  uniformly  distributed  in  depth. 

4.  The  method  consists  essentially  of  first  assuming  that  there  exists  an 
ideal  original  wave  train  at  the  locale  under  consideration  whose  spectrum 
is  obtained  using  energy  balance  equation  without  considering  wave 
breaking.  By  inposing  the  Miche  wave  breaking  criterion  (Battjes,  1974), 
an  expression  for  the  elevation  of  breaking  waves  is  established  in  terms 
of  the  elevation  of  the  ideal  original  waves  and  its  second  time 
derivative  which  are  assumed  to  be  jointly  Gaussian.  Based  on  this 
breaking  wave  model,  a  simple  but  approximate  expression  for  the  spectrum 
of  breaking  waves  is  derived  in  terms  of  that  of  the  ideal  original 
waves . 

5.  The  purpose  of  this  study  is  to  extend  the  above  method  to  the  case  in 
which  the  random  waves  are  treated  as  short-crested  whose  spectrum  has  a 
directional  spread.  The  shoreline  is  assumed  to  be  straight  and  the 
bottom  is  gently  varying  with  straight  contours  which  are  parallel  to  the 
shoreline  as  in  the  previous  study.  Two  types  of  steady,  non-uniform 
horizontal  currents  are  considered.  The  first  type  is  the  shear  current 
whose  velocity  is  parallel  to  the  shoreline  (see  Figure  2) .  The  second 
type  is  the  upwelling  type  in  which  the  current  velocity  is  normal  to  the 
shoreline.  In  each  case,  the  current  velocity  is  independent  of  the 
alongshore  coordinate  as  well  as  the  vertical  coordinate. 

6.  PART  II  describes  briefly  the  breaking  wave  model  and  gives  the 
expression  of  the  breaking  wave  spectrum  which  are  respectively 
established  and  derived  in  an  earlier  report  (Tung  and  Huang,  1967) . 

PART  III  deals  with  the  mechanics  of  wave-current  interaction  for  both 
the  shear  current  and  current  of  the  upwelling  type  without  considering 
wave  breaking.  In  PART  IV,  numerical  results  are  obtained  and  presented 
in  graphical  form.  Soam  concluding  remarks  are  given  in  PART  v. 
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PART  II.  BREAKING  NAVE  MODEL  AMD  SPECTRUM 

7 .  Ne  ASSuiM  that  wavas  braak  whanavar  tha  local  aaartical  downward 
accalaratlon  of  tha  idaal  wavaa  raachaa  a  fraction  of  tha  gravitational 
accalaration  g*  (Phillipa,  19B0) .  Mhan  this  happans*  a  portion  of  the 
mass  is  datachad  from  tha  surfaca^  tha  surfaca  elevation  is  reduced,  and 
tha  larger  the  vertical  downward  accalaration,  tha  larger  the  reduction. 

8.  Referring  to  Figure  1,  let  ^ (t)  and  ^^(t)  represent,  respectively,  the 
elevations  of  the  ideal  waves  and  the  breaking  waves  at  a  fixed  point  in 
space  where  t  is  time.  At  points  such  as  A,  where  C (t)>0  and  when 

^(t)<-Kg  (see  Equation  4  for  K) ,  wave  breaking  occurs  and  the  breaking 
wave  elevation  is  given  by 


C(t) 


(1) 


Here,  and  hereafter,  overdot  denotes  differentiation  with  respect  to 
time. 


9.  At  point  B,  where  ^(t)<0  and  when  ^(t)<-Kg,  the  breaking  wave  elevation 
is  given  by 


;b(t) 


-;(t) 

Kg 


(2) 


10.  Based  on  the  above  considerations  and  noting  that  no  wave  breaking  takes 
place  when  ^(t)>-Kg,  in  which  case  C  (t)  remains  unchanged,  C^Ct)  may  be 
written  as 


Cb  -  C  (-|^)  H(-^-Kg)H(C)+C(-^)  H(-^'-Kg)H(-C)+CH(i;+Kg)  (3) 

where  H(*)  is  the  Heaviside  unit  step  function  and,  for  brevity,  the 
argument  t  in  ^ (t)  and  ^ (t)  is  omitted. 

11.  For  linear  waves,  the  quantity  K  is  shown  in  the  earlier  study  (Tung  and 
Huang,  1987)  to  be  given  by 

K  ■  0.44tanhGijjd)  ((,) 


where  d  is  local  water  depth, 

Kq  •  w^/g  (5) 

*For  convenience,  symbols  and  abbreviations  are  listed  in  the  notation 
(Appendix  A) . 
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Is  characteristic  wave  number 
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fiii  S(u))dco 


fS  (id)  d(jJ 


(6) 


is  characteristic  frequency  and  S(to)  is  the  energy  spectrum  of  the  ideal 
random  waves  in  frequency,  ui,  space,  in  deep  water. 

12.  In  Equation  3,  the  breaking  wave  elevation  nonlinear  function  of 

2^  and  elevation  and  its  second  derivative  of  the  original  ideal 

waves  which  are  assumed  to  the  stationary  and  jointly  Gaussian  with  zero 
mean  values.  The  determination  of  the  spectrum  of  may  therefore  be 
achieved  in  a  straight  forward  manner  (Papoulis,  1965)  . 

13.  The  spectrum  of  is  obtained  by  first  forming  its  autocorrelation 
function  (t)  from  Equation  3  (where  T  is  time  lag) .  In  doing  so,  we 
have  neglected  the  second  term  in  Equation  3  based  on  the  consideration 
that  the  probability  of  occurrence  of  negative  peaks  such  as  point  B  in 
Figure  1  is  usually  small,  especially  when  the  spectrum  of  the  waves 
under  consideration  is  reasonably  narrow.  In  this  way,  the  derivation 
is  much  shortened,  and  our  conq2utation  shows  that  the  error  incurred  by 
ignoring  the  second  term  in  Equation  3  is  indeed  inq^erceptibly  small. 

14.  The  resulting  autocorrelation  function  R^j  (t)  is  a  nonlinear  function  of 

the  correlation  functions  where  Et«) 

denotes  the  expected  value  of  the  quantity  enclosed  in  the  brackets  and 
the  subscripts  1  and  2  refer  to  cjuantities  evaluated  at  time  instants 

ti  ••  t+T  and  t2“  t.  The  autocorrelation  Rj,  (l),  viewed  as  a  function  of 
these  correlation  functions,  may  be  expanded  by  Taylor's  series 
(Borgman,  1965)  .  By  retaining  only  the  zeroth  and  the  first  order  terms 
of  the  series  and  taking  the  Fourier  transform  of  the  resulting 
approximate  autocovariance  function  ,  the  approximate  spectrum  of 
breaking  waves  is  obtained  as 


S^(UJ)  -  F  (uj)S(u) 

in  which 

is  a  fourth  order  polynomial  function  of  (jDand  may  be  looked  upon  as  a 
filter  function  which  accounts  for  the  effect  of  wave  breaking  on  the 
spectrum  of  the  ideal  waves. 

15.  The  following  functions  and  parameters  that  appear  in  Equation  8  are 
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d*£ln«d  in  th«  following . 


16. 


I  _2 

Z(n)  -  exp  (  -  5  ) 


(9) 


and 


Q(C)  -  f  Z(n)dTi 

are  probability  functions  (Abramowitz  and  Stegun,1966) 

r  “  ^  cJ  duj 

(ui)  duj 


(10) 


(ID 

(12) 


and 


■  J^UJ^S  <(jj)  d(jj 

are  respectively  expected  values  E(i;^] and  E[*^2]^ 
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(13) 


(D) 


is  a  measure  of  the  extent  of  wave  breaking  which  has  been  shown  to  be 
rather  larger  than  unity  (Tung  and  Huang,  1987)  and 


1  - 


[r«>] 


r  r 


(4) 


(15) 


which  lies  between  zero  and  unity  is  known  as  the  bandwidth  parameter  of 
S(U))  (Cartwright  and  Longuet-Higgins,  1956). 


In  Equation  8, 
2 


1  > 
l»-» 

r(4) 

Ia' 

Aj  «  Bi;  +  Q(-6)  >  0 


and 


A2  *  6K  -  6Z(6)  Q  (  - 

where 

i!  ■  /  «<-  , 


B  D-e 


)- 


6Q(B/e) 

./2Tr(l-€2) 


+  SZ(6) 


n)  dn 


(16) 

(17) 

(18) 

(19) 
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PART  III:  WAVE -CURRENT-BOTTOM  INTERACTIONS 


17.  In  reference  to  Figure  2,  consider  a  straight  shoreline  and  assume  that 
the  bottom  is  gently  varying  with  contours  parallel  to  the  shore.  A 
random  short-crested  wave  system  approaches  the  shore  from  deep  water 
and  enters  a  region  of  current.  Two  types  of  horizontal  current  are 
considered.  The  first  type  is  the  shear  current  in  which  case  the 
horizontal  current  velocity  V(x)  is  parallel  to  the  shore  and 
independent  of  the  y-coordinate .  The  second  type  of  current,  referred 
to  as  the  upwelling  case,  has  its  velocity  U(z)  in  the  x-direction  and 
is  also  assumed  to  be  independent  of  the  y-coordinate.  In  each  case, 
the  current  is  steady  in  time,  slowly  varying  in  the  horizontal  space 
and  independent  of  the  vertical  coordinate. 

18.  The  wave  system  is  viewed  as  a  congregate  of  wave  conponents  with 

distinct  freq[uency  and  angle  of  approach  0.  For  each  wave  component, 
the  apparent  f requency  o) ,  in  a  stationary  frame  of  reference,  is 
invariant  and  related  to  the  intrinsic  frequency,  u3  r'  frame  of 

reference  moving  with  the  current  as 


U3  =  0)  +  U„  .  k 

r  ~H  ~ 


(20) 


where  U  is  horizontal  current  velocity  vector,  Jc  is  wave  number  vector 
and  *■“ 


u 


2 


r 


gk  tanh(kd) 


(21) 


d(x)  being  the  slowly  varying  local  water  depth,  assumed  to  be 
independent  of  the  y  -  coordinate  and  k  is  the  magnitude  of  k .  From 
irrotationality  of  wave  number  vector,  we  have 


k  sin  G  =  k  sinG 

oo  on 


(22) 


where  the  subscript  (  °°)  refers  to  quantities  evaluated  in  quiescent 
deep  water  free  of  current. 


19.  Ignoring  wave  breaking  and  using  the  energy  balance  equation  or  from 
conservation  of  wave  action  consideration,  it  was  found  (Tayfun,et  al. 
1976)  that  the  wave  frequency  spectrum  S  (uJ  ,Q )  at  the  locale  under 
consideration  is  given  by  „  . 


8(0), 0) 


k(Cc)  Ji¬ 


ll^- 

U) 


) 


where 


“.'Sr* 


(23) 


0) 


2k_ 


(24) 
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(25) 


18  the  Bagnitude  of  group  velocity  of  deep  water  waves  and 


gr 


d  u 
r 


(JL) 


is  the  magnitude  of  group  velocity  of  waves  in  relative  frame  of 
reference  and 


n 


,  2kd 

i_(l4  -  \ 

2  '  Binh  (2kd)'’ 


(26) 


for  brevity. 

20.  Thus,  for  shearing  current,  -  (  0,V(x)  ), 


uj*u  +  Vksin0 
r 


and 


S(ai,G  ) 


(i)-Vksin0 


S  (w.ej 


2k  c  +Vsin0 
“  gr 


For  the  upwelling  case,  “  (U(x),0), 


oj  =  Ci)  4-Ukcos0 

r 


and 


(27) 


(28) 


(29) 


k 

S(u),0)  »  — 


ti>“Ukcos0 


C  +UCOS0 

gr 


(30) 


21.  To  coznpute  the  spectrum  at  a  sjpecified  locale  with  given  values  of  water 
depth  d  and  current  velocity  Ujj,  and  for  a  specific  wave  component  of 
relative  frequency  and  angle  0,  the  corresponding  values  of  wave 
number  k  and  absolute  frequency  o)  are  found  from  Equations  21  and  20  (or 
Equations  27,  29)  respectively  from  which  the  wave  number  and 

angle  0<»  (from  Equation  22)  in  deep  water  are  obtained.  Finally,  S((jj,0) 
is  computed  from  Equation  23  (or  Equations  28,  30)  . 


22  . 


The  wave  components  are  all  subject  to  the  brea)cing  limit  (Tayfun  et 
al.,  1976) 


(31) 


and  the  kinematic  limit 
|sin0,|5  1 


(32) 
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Those  wave  components  which  violate  Equation  32  are  reflected  and  the 
incident  wave  field  is  accordingly  modified.  Such  modification, 
however,  is  not  carried  out  in  this  study.  Those  wave  components  which 
violate  Eqi  ^tion  31  are  cut  off  and  deleted  beyond  the  point  under 
conside'-ation.  There  also  may  occur,  as  is  the  case  in  subsequent 
numeri'^al  results,  that  for  specified  0)r,  d,  and  0>  the  energy  in 
deep  water,  is  very  s^ll  or  non-existent.  In  that  case 

S((4^  is  sot  equal  to  zero. 

23.  To  account  for  the  effect  of  wave  breaking,  the  spectrum  S(u),0),  viewed 
as  a  frequency  spectrum  for  specified  value  of  0,  is  first  transformed 
into  the  frame  of  reference  moving  with  current.  That  is 


S((ij  ,0)  =  S((jo,0) 

T  du 

r 


(33) 


where  du/dwr  is  obtained  from  Equation  20  (or  Equations  27,  29).  S(ui^,Q) 
is  then  used  in  place  of  S  (oj)  in  Equation  7  from  which  the  breaking  wave 
spectrum  in  relative  frame  of  reference  ,0)  is  determined.  Spectrum 

of  breaking  waves,  Sjj(u),0),  in  stationary  frame  of  reference,  is  then 
calculated  by  transformation.  That  is. 


S^(u.0) 


du) 


(3A) 


24.  Although  the  spectra  of  the  breaking  waves,  in  relative  and  absolute 
fraunes  of  references  so  obtained  as  outlined  above  are,  within  the 
confines  of  the  assmqstions  made,  strictly  correct,  iit^srovements  of  the 
results  are  nevertheless  possible  by  carrying  out  the  computations 
iteratively.  That  is,  upon  obtaining  the  breaking  wave  spectrum 

in  the  relative  frame  of  references,  it  may  be  used  as  the  ideal 
original  wave  spectrum  S (m)  in  Equation  7  from  which  a  revised  spectrum 
Sb(‘^r'"^  results.  This  process  may  be  repeated  until,  say,  the  peak 
value  of  converges  to  within  a  prescribed  limit  of  tolerance. 

The  breaking  wave  spectrum  S^tJ,0)  in  the  stationary  frame  of  reference 
is  then  computed  according  to  Equation  34.  In  this  study,  the  tolerance 
limit  is  set  at  0.001  m^/s  for  all  cases  considered.  The  maximum  number 
of  cycles  of  iteration  necessary  to  satisfy  this  tolerance  limit  is 
eight . 

25.  The  iteration  process  described  above,  although  intuitively  appealing 

and  presents  no  problem  in  confutation,  is  not  without  difficulty  in 
theoretical  terms  since  Equation  7  is  derived  based  on  the  assumption 
that  the  ideal  waves  are  a  zero  mean  Gaussian  process.  However,  in  the 
earlier  study  (Tung  and  Huang,  1987),  it  was  found  that  the  mean  value 
of  is  rather  small  and  in  another  study  (Tung  et  al.,  1966)  the 
deviation  of  the  probability  function  of  C b  Gaussian  is  also 

minimal  so  that  the  iterative  scheme  adopted  in  this  study  must  be 
considered  acceptable. 
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PART  XV:  NUMERICAL  RESULTS 

26.  In  order  to  obtain  numerical  results,  the  directional  wave  spectrum  in 
deep  water  mu^kt.  be  specified.  In  this  study,  it  is  assumed  that  the 
directional  wave  spectrum  takes  the  form 


S(uj.ej  -  S^(co)t(GJ 


(35) 


where 


3Tr 


cos 


(36) 


is  the  directional  energy  spreading  function.  For  5^(0)),  there  are  a 
great  number  of  frequency  spectra  (The  JONSWAP  Spectrum,  for  example) 
one  can  choose  from.  In  this  study,  we  use  the  wallops  spectrum  (Huang 
et  al,  1981)  whi  .h  takes  the  form 


27  . 


s«(^) 


2 

og _ 

m  5-m 

u 


(37) 


The  quantity  m  gives  the  magnitude  of  the  slope  of  the  frequency 
spectrum  (on  log-log  scale)  in  high  frequency  range  and  is 


l0g(2TT^§S 

log2 


(38) 


where 


(39) 


1 

o 

is  the  significant  slope  of  the  wave,  \  •2'>T/k  being  the 
characteristic  wave  length  (see  Equatioi9  5) .  °rhe  quantity  Q  is  given  by 


a 

where 


(ffl-D/A  , 

^ _  (2710^ 

“  ^(m-5)/4  r(  (m-l)/4) 

r  {•)  is  the  gamma  function 


(Abramowitr  and  Steg 


(40) 
1968)  . 


28.  The  Wallops  spectrum  is  thus  seen  to  depend  on  two  .T>eters,  the 

significant  slope  §  ,  which  defines  m  and  hence  d,  -  -  o'  ‘'^^ch  is  the 
frequency  corresponding  to  the  peak  of  the  Wallops  spectr^.  m  what 
follows  all  numerical  results  are  obtained  for  values  of  0.015  and 

-  0.6  rad./s. 

29.  The  directional  spectrum  S  in  deep  water  as  given  by  Equation  34 

is  shown  graphically  in  two  ways.  In  Figure  3,S^  are  given  for 

various  values  of  0^  as  function  of  w.  Since  S  (‘*',^35).^  is  symmetrical 
about  0  ■  0,  the  spectra  are  shown  only  for  0  <  0^^  <7  •  in  Figure  4, 
the  contour  lines  of  S„(  w,0^)  are  dtawn.  The  ordinates  of  the  contour 
lines  are  marked  as  sho%m. 
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30.  In  the  following,  we  first  consider  the  case  in  which  deep  water  wave 

system  with  S^(u),6^)  shown  in  Figures  3  and  4  meets  a  shear  current  of  v 
-  2  deep  water.  Qualitative  behavior  of  wave  con^ionents  on  a 

positive  shear  current  in  deep  water  is  shown  in  Figure  5a  and  5b.  The 
ideal  wave  spectr\un  S(u, 6)  and  the  breaking  wave  spectrum  S]^(u>,0)  are 
given  in  absolute  frame  of  reference  for  discrete  values  of  G  in  Figures 
6  to  14.  The  contour  lines  of  S^{u),Q)  in  this  case  are  shown  in  Figure 
15.  From  Figures  6  to  14,  it  is  seen  that  most  of  the  wave  energy  is 
concentrated  near  small  values  of  0  and  no  wave  breaking  occurs.  For 
waves  with  0  -  -40°,  and  -30°,  the  spectra  terminate  at  frequencies 
1.205  and  2.208  rad/s  respectively.  Naves  of  frequencies  beyond  these 
cut-off  frequencies  originate  from  regions  of  (u  spsce  where  no 

wave  energy  exists  (see  Figure  5b) .  Naves  with  0<O  are  subject  to 
limits  set  by  wave  breaking  but  the  breaking  limits  are  much  larger  than 
the  cut-off  frequencies  shown  above.  Likewise,  waves  with  Oo'are 
subject  to  kinematic  limits  and  may  be  reflected.  However,  these 
kinematic  limits  are  well  beyond  frequencies  where  energy  of  component 
waves  of  any  significance  exists. 

31.  We  next  consider  the  case  in  which  deep  water  waves  of  spectra  shown  in 

Figures  3  and  4  encounter  an  adverse  current  of  the  upwelling  type  with 
U  — 2"'/g  while  in  deep  water.  Qualitative  behavior  of  waves  on  adverse 
current  in  deep  water  is  sketched  in  Figure  16.  The  ideal  wave  spectrum 
S((jj,G)  and  the  breaking  wave  spectrum  0)  in  absolute  frame  of 

reference  are  given  in  Figures  17  to  21,  for  various  values  of  0.  The 
contour  lines  of  St)(ai0)  are  given  in  Figure  22.  Because  of  synsnetry, 
only  spectra  and  contour  lines  for  G>0®  are  shown.  As  may  be  expected, 
wave  energy  is  highest  at  0-  0°.  For  waves  with  0-0°  and  10°,  wave 
breaking  limits  the  frequencies  at  w-  1.225  and  1.244  rad/s 
respectively  where  wave  energy  grows  indefinitely.  In  reality,  waves 
break  long  before  these  limits  are  reached.  For  waves  with  0-  20°, 

30°,  and  40°,  waves  beyond  u-  1.262,  1.16,  and  1.00  respectively  do  not 
exist  since  these  originate  from  regions  of  (u)-G^)  space  where  there  is 
no  wave  (see  Figure  16) .  For  0-0°,  10°  and  20°,  wave  breaking  reduces 
the  ideal  wave  spectrum  S((4Q),  but  beyond  0=  30°,  no  wave  breaking  is 
seen  to  occur.  Other  than  for  the  cases  of  0-0°  and  10°,  the  limits 
set  by  wave  breaking  are  beyond  those  cut-off  frequencies  shown  in  these 
Figures.  From  Figure  22,  it  is  seen  that  the  spectrum  Sb(‘^,  0)  is  double 
peaked. 

32.  When  deep  water  waves  with  spectrum  given  in  Figures  3  and  4  propagate 
toward  the  shore  in  still  water  without  current  over  a  gently  varying 
sea  bed  with  straight  contours  parallel  to  the  straight  shoreline,  the 
waves  are  re^ -acted  and  bend  toward  the  shore  as  shown  in  Figure  16.  At 
water  depth  -  10*",  the  spectra  S(a),G)  and  Sb(uvG)  are  computed  and 
given  for  rete  values  of  ©in  Figures  23  to  27  and  the  contour  lines 
of  Sb((jj,  6)  are  given  in  Figure  28.  Due  to  symmetry,  only  cases  of  G>0°, 
are  given.  M'''"  of  the  wave  energy  is  contained  in  the  neighborhood  of  0 
-  0°.  F(  T  0 »  0°,  10°  and  20°,  wave  breaking  reduces  wave  energy,  but 
for  waves  with  0-  30°  and  40°,  the  effect  of  wave  breaking  is 
vanishingly  small.  For  0-  10°,  20°,  30°,  and  40°  there  exist  cut-off 
frequencies  u)  -  0.197,  0.370,  0.521,  and  0.695  respectively  below  which 
no  wave  exists.  These  low  frequency  long  waves  which  are  moat  severely 
affected  by  the  bottom,  originate  from  regions  of  (a)-0^)  space  where 
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ther«  is  no  «av«.  High  fr«qu«ncy  short  waves  are  essentially  deep  water 
waves  iu)af£ected  by  the  bottom. 

33.  We  now  consider  the  case  in  which  the  deep  water  wave  system  with 
spectrum  shown  in  Figures  3  and  4  enters  a  region  of  shear  current  of 
velocity  V  -  2^/ g  approaching  a  straight  shore  over  gently  varying 
seabed  of  straight  contours  parallel  to  the  beach.  At  water  depth  d  ■> 

10  the  spectra  SCjJ,6)  and  S)3<bv6)  in  absolute  frame  of  reference  are 
given  in  Figures  29  to  35  and  the  contour  lines  of  S^((i),Q)  are  given  in 
Figure  36.  Similar  to  the  cases  of  waves  on  shear  current  in  deep  water 
and  waves  approaching  the  shore  without  current,  wave  energy  is  highest 
around  0  -  0°.  Nave  breaking  occurs  for  waves  with  0  -  -20°,  -10°,  0°, 
and  10°,  but  not  for  waves  with  larger  values  of  G.  Figure  37  gives  a 
sketch  of  behavior  of  waves  on  positive  shear  current  in  water  of  finite 
depth  with  parallel  contours.  Long  waves  are  affected  predominantly  by 
the  bottom  whereas  short  waves  are  affected  by  the  current.  For  waves 
with  0<O°,  both  short  and  long  waves  are  bent  toward  the  shore,  but  for 
waves  with  Q>0°,  short  waves  bent  by  the  current  become  increasingly 
parallel  to  the  shore  but  refraction  by  the  bottom  bends  the  long  waves 
toward  the  shore.  Thus,  in  Figures  29  to  35,  for  both  waves  with  G<0° 
and0  >0°,  there  are  lower  cut-off  frequencies  of  <^“0.609,  0.401,  0.190, 
0.197,  0,204,  0.364,  and  0.477  rad/s  corresponding  respectively  to 

cases  of  0  -=  -30°,  -  20°,  -10°,  0°,  10°,  20°,  and  30°.  Naves  of 
frequencies  below  these  cut-off  frequencies  originate  from  regions  of 
(u)-0o4  space  where  no  wave  is  present.  The  short  waves  travel 
essentially  in  deep  water  on  the  shear  current  and  behave  much  the  seme 
as  those  given  in  section  3  where,  for  waves  with  0  “  -30°,  there  is  an 
upper  cut-off  frequency  of  u*  2.208  rad/s. 

34.  The  last  case  considered  is  that  of  deep  water  waves  of  spectra  shown  in 
Figures  3  and  4  encountering  a  current  of  the  upwelling  type  with  U  *  - 
2"'/g  in  finite  water  depth.  For  waves  at  water  depth  d  «  lO’'",  the 
spectra  S(u),G)  and  (av0)  in  absolute  frame  of  reference  are  given  in 
Figures  38  to  41  and  the  contour  lines  of  Sjj((u,G)  are  given  in  Figure 
42-  Because  of  symmetry,  only  the  cases  of  Q>0°  are  given.  Energy  is 
concentrated  around  0**  0°  and  wave  breaking  occurs  for  waves  with  0« 
0°,  10°,  and  20°,  but  not  for  waves  with  0  -  30°  much  the  same  as  in  the 
case  of  deep  water  waves  meeting  an  upwelling  current  in  deep  water. 

Both  short  and  long  waves  are  bent  toward  the  shore  with  the  short  waves 
affected -by  the  current  and  the  long  waves  by  the  sea  bottom.  Thus,  for 
the  short  waves,  there  are  upper  cut-off  f recjuencies,  as  limited  by  wave 
breaking  of  oj “  1.225  and  1.244  rad/s  corresponding  to  waves  with  G  “  0° 
and  10°  respectively.  For  waves  with  0  ■  20°  and  30°,  the  upper  cut-off 
frequencies  are  respectively  u  -  1.262  and  1.15  rad/s,  as  in  the  case  of 
waves  encountering  an  upwelling  current  in  deep  water.  For  the  long 
waves,  the  lower  cut-off  frequencies  are  U)  ■  0.234,  0.450  and  0.739 
rad/s  corresponding  respectively  to  the  cases  of  0  »  10°,  20°  and  30°. 
Waves  of  frequencies  below  these  values  originate  from  regions  of  (u-Goo) 
space  where  no  wave  exists. 
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PART  V:  CONCLUSION 


35.  In  an  earlier  report  (Tung  and  Huang,  1987)  we  introduced  a  method  by 
means  of  which  the  effect  of  wave  breaking  or  wave  spectrum  can  be 
incorporated.  The  method  was  applied  to  the  ease  in  which  a  uni¬ 
directional  deep  water  wave  ayatem  propagates  toward  the  shorw  over 
gently  varying  sea  bed  of  straight  parallel  contours  encountering  an 
adverse  current,  in  this  report,  the  method  is  extended  to  the  case  in 
which  the  deep  water  wave  system  represented  by  a  directional  spectrum 
propagating  over  a  sea  bed  having  the  aasm  characteristics  meeting 
current,  of  either  the  shearing  type  parallel  to  the  shore  or  of  the 
upwelling  type  perpendicular  to  the  shore. 

36.  The  method  of  incorporating  wave  breaking  is  approximate,  but  easy  to 
apply.  The  solution  to  the  wave-current  interaction  problem  (without 
considering  wave  breaking)  is  for  the  specific  case  where  the  sea  bed 
has  parallel  contours.  The  results  so  obtained  are  therefore  restricted 
in  application.  The  method  of  dealing  with  wave  breaking  can  be 
improved,  but  such  refinement  is  not  considered  necessary  at  this  stage 
of  development.  The  restriction  imposed  by  the  idealization  of  the  sea 
bed  can  be  relaxed  by  adopting  a  fully  numerical  treatment  of  the  wave- 
current  interaction  problem  (without  considering  wave  breaking) . 

37 .  The  results  of  this  study  must  be  checked  against  those  of  other 
analytical  methods  and  experimental  data,  both  from  the  field  and 
laboratory. 
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Figure  1.  Wave  profile 


Figure  2.  Definition  sketch  of  wave,  current 
and  bottom  contour 
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Figure  3.  Spectra  Soo(w,G<x,)  of  deep  water  waves  under 

zero  current  condition,  for  0.^  =  0°,  lO",  20°,  30°,  40°, 

and  50°,  and  for  §  =  0.015  and  w  =  0.6  rad/sec 
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Figure  4.  Contour  lines  of  S„(u),0)  of  deep  water 

waves  under  zero  current  condition,  for  §  =  0.015 

and  0)  =0.6  rad/sec 
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Figure  5. 
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Qualitative  behavior  of  deep  water  waves 
on  positive  shear  current 
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Figure  6.  Spectra  S(a),G)  and  Sjj(w,0)  of  deep  water 
waves  on  current  with  V  =  2™/g  ,  for  0  =  -40°,  and  for 
§  =  0.015  and  u)  =0.6  rad/sec 
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Figure  9.  Spectra  S(w,0)  and  SwCto.Q)  of  deep  water 
waves  on  current  with  V  =  2™/g  ,  for  0  =  -10°,  and  for 
§  =  0.015  and  o)  =  0.6  rad/sec 


0  0.2  0.4  0.6  0.3  :  1.2  1.4  1.6  1.8  I 

Figure  10.  Spectra  S(a},0)  and  S^{ni,Q)  of  deep  water 
waves  on  current  with  V  =  2”^/g  ,  for  0  =  0°,  and  for 
§  -  0.015  and  w  =0.6  rad/sec 
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Figure  11.  Spectra  S(tu,0)  and  S^^Coj,©)  of  deep  water 

waves  on  current  with  V  =  2"*/s  .  for  0  =  10°,  and  for 

§  =  0.015  and  o)  =  0.6  rad/sec 
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Figure  13.  Spectra  S(u,0)  and  Sb(w,0)  of  deep  water 
waves  on  current  with  V  =  2™/s  ,  for  0  =  30°,  and  for 
§  =  0.015  and  w  =0.6  rad/sec 
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Figure  14.  Spectra  S(a),0)  and  S|j(cj,0)  of  deep  water 
waves  on  current  with  V  =  2®/  ,  for  0  =  40°,  and  for 

§  =  0.015  and  oi  ^=  0.6  rad/sec 
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Figure  16.  Qualitative  behavior  of  deep  water 
waves  on  adverse  current  of  upwelling  type  and 
of  waves  over  bottom  of  parallel  contours 
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Figure  17.  Spectra  S(u),0)  and  Sjj(u),0)  of  deep  water 

waves  on  current  with  U  =  -2"'/g  >  for  0  =  0°,  and  for 

§  =  0.015  and  u  -  O.C  rad/sec 
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Figure  18.  Spectra  S(u),0)  and  Sjj(u),0)  of  deep  water 
waves  on  current  with  U  =  -2®/g  ,  for  0  =  10  ,  and  for 
§  =  0.015  and  o)  =  0.6  rad/sec 


Figure  19.  Spectra  S(a),0)  and 
waves  on  current  with  U  =  -2*"/g 
§  =  0.015  and  w 


Sj3(a),0)  of  deep  water 
,  for  0  =  20° ,  and  for 
=  0.6  rad/ sec 
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Figure  20.  Spectra  S(OJ,0)  and  Sjj((xj,0)  of  deep  water 

waves  on  current  with  U  =  -2®/g  ,  for  0  =  30°,  and  for 

§  =  0.015  and  w  =  0.6  rad/sec 
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Figure  21.  Spectra  S(u),0)  and  Sjj(u),0)  of  deep  water 
waves  on  current  with  U  =  “2^/g  ,  for  0  =  40° ,  and  for 
§  =  0.015  and  w  =  0.6  rad/sec 


Figure  22.  Contour  lines  of  S{j(to  0)  of  deep 

water  waves  on  current  with  U  =  -2"'/g  ,  for 

§  =  0.015  and  u  =  0,6  rad/sec 
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Figure  23.  Spectra  S(oJ,0)  and  S}j(a),0)  of  waves  at  water 
depth  d  =  10™  under  zero  current  condition,  for  0=0®, 
and  §  =  0.015  and  u  =0.6  rad/sec 
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Figure  24.  Spectra  S(w,0)  and  Sjj(w,0)  of  waves  at  water 
depth  d  =  10™  under  zero  current  condition,  for  0  =  10®, 
and  §  =  0.015  and  w  =0.6  rad/sec 
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Figure  25.  Spectra  S(a),0)  and  S^,(u),0)  of  waves  at  water 
depth  d  =  10™  under  zero  current  condition,  for  0  =  20°, 
and  §  =  0.015  and  w  =0.6  rad /sec 
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Figure  26.  Spectra  S((jj,0)  and  S(j(w,0)  of  waves  at  water 
depth  d  =  10™  under  zero  current  condition,  for  0  =  30°, 
and  §  =  0.015  and  o)  =  0.6  rad/sec 
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Figure  27.  Spectra  S((o,0)  and  S^j(to,0)  of  waves  at  water 
depth  d  =  10™  under  zero  current  condition,  for  0  =  40“, 
and  §  =  0.015  and  o)^  =  0.6  rad/sec 


Figure  28.  Contour  lines  of  Sij(u),0)  of  waves  at 
water  depth  d  =  10™  under  zero  current  condition, 
for  §  =  0.015  and  u)  “0.6  rad/sec 
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Figure  34.  Spectra  S(a),0)  and  S{j((i),G)  of  waves  at  water 
depth  d  =  10"'  on  current  with  V  =  2"'/g  ,  for  0  =  20°,  and 

for  §  =  0.015  and  0)  =0.6  rad/sec 
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Figure  35.  Spectra  S(oj,0)  and  for  waves  at  water 

depth  d  =  lO™  on  current  with  V  =  2"'/g  ,  for  0  =  30° ,  and 
for  §  =  0.015  and  u)  =0.6  rad/sec 


40 


c/3 

B 


Figure  40. 
depth  d  = 


Crt 

«N 

B 


Figure  41. 
depth  d  = 


Spectra  S(a),0)  and  S^j(w,0)  of  waves  at  water 
10®  on  current  with  U  =  -2®/g  ,  for  0  =  20°,  and 
for  §  =  0.015  and  w  =  0.6  rad/sec 


Spectra  S(w,0)  and  S{,(u),0)  of  waves  at  water 
10®  on  current  with  U  *  ~2®/g  ,  for  0  =  30°,  and 
for  §  =  0.015  and  U)  =  0.6  rad/sec 


Figure  42.  Contour  lines  of  S^(a),G)  of  waves  at 

water  depth  d  =  lO*"  on  current  with  U  =  -2"’/g  , 

for  §  =  0.015  and  u  =  0.6  rad/sec 
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Quantities  defined  in  Equations  17  and  It,  respectively 
Group  velocities  defined  in  Equation  24  and  25,  respectively 
Water  depth 

Expected  value  of  the  quantity  enclosed  in  the  brackets 

Filter  function  defined  in  Equation  8 

Gravitational  acceleration 

Heaviside  unit  step  function 

Coefficient  defined  in  Equation  4 

Magnitude  of  wave  nuoiber  vector  k 

Nave  niinber  vector 

Magnitude  of  wave  number  in  deep  water 
Characteristic  wave  number  in  deep  water 

Magnitude  of  the  slope  of  the  Nallops  spectrum  for  deep  water 
waves  in  the  high  frequency  range  on  log-log  scale  given  in 
Equation  38 

Quantity  defined  in  Eq[uation  19 
Quantity  defined  in  Equation  26 
Probability  function  defined  in  Equation  10 
Autocorrelation  function  of  breaking  wave  elevation 
Quantities  defined  in  Equations  11,  12,  and  13,  respectively 

Ideal  wave  spectrum 
Breaking  wave  spectrum 
Spectrum  of  waves  in  deep  water 

Directional  spectrum  of  ideal  waves  in  sbsolute  fraas  of 
reference 

Directional  spectrum  of  breaking  waves  in  absolute  fraxae  of 
reference 
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Directional  spectrua  of  ideal  waves  in  relative  frame  of 
reference 

Directional  spectrum  of  breaking  waves  in  relative  frame  of 
reference 

Directional  Spectrum  in  quiescent  deep  water 
Time 

Time  instants  t  +  T  and  t»  respectively 

Speed  of  current  of  the  upwelling  type  in  the  x-direction 

Horizontal  current  velocity  vector  and  its  magnitude, 
respectively 

Speed  of  shear  current  in  the  y-direction 
Horizontal  coordinates  defined  in  Figure  1 
Probability  function  defined  in  Equation  9 
Coefficient  defined  in  Equation  40 
Wave  Breaking  parameter  defined  in  Equation  14 
The  gamma  function 

Spectral  bandwidth  parameter  defined  in  Equation  15 
Elevations  of  ideal  and  breaking  waves,  respectively 
Dummy  variable 

Angles  of  wave  at  locale  under  consideration  and  in  deep  water, 
respectively 

Characteristic  wave  length 
Dummy  variable 


T  Time  lag 

$(•)  Directional  energy  spreading  function  defined  in  Equation  36 

b)  Wave  frequency 

Wave  frequency  In  relative  frame  of  reference 

b>  Characteristic  wave  frequency  defined  in  Equation  6 
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(Dq  P«raMt«r  of  Wallops  Spsetrun 

(i>2^  Quantity  dsflnad  in  Equation  IS 

Subactiots 


1,2 


Quant itias  avaluatad  at  tis»  instants  tx«  snd  t2,  raspaetivaly 
Quant it las  avaluatad  in  daap  watar  undar  saro  currant  condition 


Symbols 


§ 


Diffarantiation  with  raspact  to  tina 
Significant  wave  slopa  dafinad  in  Equation  39 
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